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ABSTRACT 

Recent far-infrared (IR) observations of supernova remnants (SNRs) have revealed significantly large 
amounts of newly-condensed dust in their ejecta, comparable to the total mass of available refractory 
elements. The dust masses derived from these observations assume that all the grains of a given species 
radiate at the same temperature, regardless of the dust heating mechanism or grain radius. In this 
paper, we derive the dust mass in the ejecta of the Crab Nebula, using a physical model for the heating 
and radiation from the dust. We adopt a power-law distribution of grain sizes and two different dust 
compositions, and calculate the heating rate of each dust grain by the radiation from the pulsar wind 
nebula (PWN). We find that the grains attain a continuous range of temperatures, depending on their 
size and composition. The best-fit model to the observed IR spectrum consist of amorphous carbon 
grains with a total mass of 0.027 ± 0.003 M0. We find that the power-law size distribution of dust 
grains is characterized by a power-law index of 3.5 and a maximum grain size larger than 0.6 /im. 
The grain sizes and composition are consistent with what is expected for dust grains formed in a 
Type IIP SN. Our derived dust mass is significantly less than the 0.11-0.24 Mq reported in previous 
studies of the Crab Nebula that are based on a simplified two-temperature models. These models 
also require a larger mass of refractory elements to be locked up in dust than was likely available in 
the ejecta. The results of this study shows that a physical model resulting in a realistic distribution 
of dust temperatures can significantly affect derived dust masses. Our study has therefore important 
implications for dust mass estimates in other SNRs and for the ultimate question of whether SNe are 
major sources of dust in the Galactic interstellar medium (ISM) an in external galaxies. 
Subject headings: dust, extinction - infrared: ISM - ISM: individual objects (Crab Nebula) - ISM: 
supernova remnants - pulsars: individual (PSR B0531-I-21) 
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1. INTRODUCTION 

In recent years, the far-infrared (IR) observations of 
supernova remnants (SNRs) with the Herschel Space Ob- 
servatory have allowed us to finally tackle the longstand- 
ing question of whether supernovae (SNe) contribute a 
significant amount of dust to the interstellar medium 
(ISM). Theoretical dust condensation models predict 
that 0.1-0.7 MfD of dust should form in a SN explo- 
sion f Todini fc Ferrara"200r; 'Cherc hneff fc Dwekl ^OWt 
iNozaw a ct al. 2010; Silvia ct al. 20 ij), with ~ 40% sur- 
viving the reverse sho ck and beiiig injected into the 
surro unding ISM (e.g., iDwek et"an 120081 : iKozasa et al.l 
120091 ). Observations of dusty galaxies at high redshifts 
also imply that large quantities of dust are produced on 
short timescales, requiring 0.1-1.0 M0 of dust to be pro- 
duced by an average SN e xplosion, dep ending on the dust 
lifetime in the ISM fe.g.. iDwek et al.1.2009' ). 

Spitzer Space Telescope observations of SNRs have 
revealed ejecta dust in many SNRs with estimated 
masses in the 0.02-0. IMrr^ range (iSugermaii et al. 120061 : 
IRho et all [2008L I200I iTemim et al.l I20ia I2012bl) . but 
the recent far-IR observations that are sensitive to cooler 
dust are now suggesting even larger masses. Gas A ap- 
pears to contain 0.075 Mq of cool (~ 35 K) ejecta dust 
located on interior of the reverse shock (jBarlow et al.l 
120101 iSibthorpe enil 120101 ). Based on recent Herschel 
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observations. lMatsuura et al.l (120111) reported 0.4-0. 7 M© 
of cool dust in SN 1987A, while iGomez et al.l (|2012[ ) find 
as much as 0.24 Mq of cool dust in the Grab Nebula. 

While these recent results imply that SNe may indeed 
be major suppliers of dust, the masses in each each of 
these cases are derived using somewhat simplified models 
that consist of one or two dust temperature components. 
A more realistic scenario is to identify the heating source 
in the ejecta, and use a continuous size distribution of 
grains with each grain heated to a different temperature 
that depends on its composition, size and optical prop- 
erties. 

In this work, we present a detailed model of the phys- 
ical mechanism giving rise to the observed IR emission 
from the Grab Nebula. Dust in the Grab Nebula was dis- 
covered as an IR excess above the synchrotron power-law 
spectrum of its pulsar wind nebula (P WN) (Trimble 1973 
Glaccum ct al. 1982; Marsden et al. 1984; Dou vion et all 
2001: Green ctal. 2004; Tcmim et al. 2006). Absorp- 
tion features from dust arc o bserved to spatially coincide 



with the ejecta fila.ments (.Woltier fc Veron-Gettv 
Fesen fc Bla"hlll990l : iHester et al.lllQQOl: iBlair et al 



1987 



1997 



LolKOlO"). and a recent analysis of the Spitzer IRS spec 



tra confirm ed that the dust is in deed located in the fila- 
ment cores (jTemim et al.ll2012b[ ) , and therefore dust that 
condensed in the ejecta. 

The model considers three distinct dust compositions, 
each characterized by a power-law distribution in grain 
radii. The primary goal of our study is to determine 
whether a physical model for the IR emission, character- 
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Fig. 1. — Absorption coefficients as a function of wavelength 
for siU cate and amorphous carbon (BE) grains from Zub ko et al.l 
l|2004l ) are shown in red and g r een, and amorphous carbon (AC) 
grains from lRouleau &: Martini 1)19911 ) in blue. 

ized by a continuous distribution in dust temperatures 
will affect the derived dust mass in this remnant. Our 
choice of the Crab Nebula is primarily motivated by the 
fact that the dust grains are heated by the synchrotron 
radiation from the PWN, so that their temperature can 
be accurately derived for each size and composition. Fur- 
thermore, since the newly-formed dust has neither been 
processed by the reverse shock, nor mixed with the ambi- 
ent medium, the Crab Nebula offers a unique opportunity 
to study the mass, composition, and size distribution of 
pristine SN-condensed dust. 

In spite of the fact that we concentrate on modeling 
the IR emission from the Crab Nebula, our results have 
global implications for deriving dust masses from IR ob- 
servations, namely, that single or two (warm, cold) tem- 
perature representations of spectra can results in a sig- 
nificant overestimate of dust mass. 

2. DUST HEATING MODEL 

2.1. Heating by the PWN 

Since the dust in the Crab Nebula is concentrated in 
the Ray lei gh- Taylor filaments that form a cage around 
the PWN (jTemim et al.ll2012b[) . we model the heating of 
dust by assuming that the heating source is located at the 
center of th e Crab Nebu l a, at th e location of the Crab 
pulsar. In iTemim et al.l (|2012b( ). we showed that the 
dominant heating source for the dust in the Crab Neb- 
ula is the synchrotron radiation from the PWN, with an 
insignificant possible contribution from collisional heat- 
ing by the gas in the filaments. In this case, the heating 
rate of a single dust grain of radius a is given by 

Hi-) = , (1) 

where is the non-thermal specific luminosity of the 
Crab Nebula's PWN, a is the grain size, Q{v,a) is the 
absorption coefficient for a given grain composition, and 
r is the distance between the radiation source and the 
dust grain. The non-thermal luminosity Li,{v) of the 
Crab Nebula's PWN th at was used in our model is sum- 
marized in F i gure 2 of iHesteii ()2008D . As described in 
ITemim et all (|2012b[) . we consider Lv{v) up to an en- 
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Fig. 2. — Best-fit dust grain heating models for three different 
grain compositions (see Figure [T]l. The da ta includes the aver- 
age Spitzer IRS dust spectrum I jTem im et al. 2012b), and the dust 
fluxes from Herschel PACS and SPIRE (Gomez et al. 2012). The 
best-fit parameters are summarized in Tabic [l] and corresponding 
total dust masses in Table 

ergy of about 0.6 keV, since we find that the fraction 
of the energy deposited in the dust at higher energies 
makes a negligible contribution of less than 1 %, due to 
the the combined effects of decreasing of both L^{v) with 
energy, and the efficiency of the energy de posited by th e 
photoelectrons in the dust grain (Dwck fc Smith|[T996[ ). 

Based on the three-dimensional models of the Crab 
Nebula, the ejecta filaments are located between 0.55-1.0 
pc from the center of the nebula (e.g. lCadez et al.ir2004l) . 
In our dust model, we allow the distance from the heating 
source r to vary between 0.5-1.5 pc in intervals of 0.2 
pc. We allow the upper limit on the distance to extend 
beyond the physical location of the filaments in order to 
account for any attenuation by dust that would affect the 
heating rate. The best-fit models described in Section [3] 
favor the low end of the distance range, suggesting that 
the internal a bsorpt ion in the nebula is negligible. In 
ITemim et all (|2012bD . we find that the optical depth is 
indeed low, with r < 1. 

In equilibrium, the heating rate in Equation ([1]) of each 
grain of radius a is equal to the radiative cooling rate 
given by 

Lgr{a) = 4™^ j TTB^iy, T)Q{v, a)dv, (2) 

where is t he Planck function and T is the dus t grain 
temperature (|Dwek fc Werneii 119811: iDwek et al.ll2OO80 . 
We used this relationship to compute the temperature 
for each grain radius and composition, and each distance 
r, and produced a final set of dust emission models by 
convolving the resulting emission spectrum for each grain 
size with the size distributions functions described in the 
previous section. 

2.2. Dust Composition 

The IR sp ectrum of dust in th e Crab Nebula is fairly 
featureless (jTemim et al.l l2012b( ). consistent with the 
generally featureless spectra generated by silicate and 
carbon grains. Theoretical dust condensation models do 
indeed show that the largest fraction of dust formed in 
Type Hp SNe is in the form of silicate and carbon (e.g. 
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Fig. 3. — The range of sp ectra simulated by the Monte Carlo 
method described in Section 12.51 is shown as the gray band, with 
the actual data overlaid in black. 

iKozasa et a"l1 [2009:l. 

We use three different sets of optical constants to cal- 
ctilate the absorptions coeffic ients Qabs{v,a): those from 
IWeingartner fc Dra inc (2001) to characterize the sihcate 
grains; those from Rouleau & Martin (1991) to charac- 
terize the amorphous carbon grai ns (labeled AC thro ugh- 
out the paper); and those from iZubko et al.| ()2004i) for 
amorphous carbon grains of type BE. The absorption co- 
efficient as a function of wavelength for each of the three 
grain compositions is shown in Figure m We note that 
Qabs{v: a) for the BE carbon does not have coverage be- 
low ~ 0.1 ^m, where a significant fraction of the energy 
absorption takes place. We only included them in this 
work fo r co mparison with previo us work bv lTemim et alj 
(|2012b( ) and lGomez et al.l(|2012l ). We later show that the 
lack of short wavelength coverage for the optical con- 
stants significantly impacts the results of the dust heat- 
ing model. 

2.3. Grain Size Distributions 

IR observations of unmixed young remnants, such as 
Gas A or SNR 1987A, have been used to derive the dust 
composition and mass, since determination of the grain 
size distribution requires detailed knowledge of the phys- 
ical condition in the ejecta, and the dust heating mech- 
anism. The size distribution of SN condensed dust is 
therefore a priori unknown. Models for the formation o f 
dust in SNe (|Todini fc Ferraral[200TI : fNozawa et al.ll2010( l 
have derived the size distribution of the various dust 
species that formed in the ejecta. However, these cal- 
culations assumed a uniform ejecta, and that the growth 
of the dust grains proceeded only through the accretion 
of single monomers, leaving out possible growth through 
coagulation in the ejecta. 

onsidering the above mentioned uncertainties, we 
adopt a general parametrization for the grain size dis- 
tribution, described by a power law in grain radii, a~", 
and a lower and upper grain radius cutoff on the grain 
radii, Omm and Umax resp ectively. His t orical ly, such 
characterization was used bv lMathis et"all (|l97l (MRN) 
to derive the size distribution of interstellar dust from 
the observed UV-o ptical extinction (for a review see 
iGlavton et al.l [2"003[ ). The incorporation of additional 
observational constraints: the diffuse IR emission, in- 



FlG. 4. — Grain temperature as a function of grain size for dust 
in the Crab Nebula. The input spectrum of the heating source 
is the broadband non-thermal spect rum of the Crab Nebula, as 
summarized in Figure 2 of HestiB 1(2001). The dust grains are 
found to be located 0.5 pc from the center of the PWN. The colored 
curves correspond to grains with absorption coefficients shown in 
Figure [T] 

terstellar abundances, and diffuse interstellar scattering, 
have yie lded a more complex interste l lar grain size distri- 
bution (IWeingartner fc Draind [200lt iLi fc Draind 120011: 
IZubko et al.ll2004D . This distribution is determined by 
the size distribution of the grains that are injected from 
all sources into the ISM, and by the various interstel- 
lar processes that alter their sizes, including thermo- 
kinetic sputtering, evaporative and shattering grain- 
grain collisions, accretion and coagulation. Attempts to 
characterize the net size distribution res ulting from all 
these processes were made by jLi ffman fc Clavtoiilll989l : 
lO'Donnell fc MathisllT997l: IHirashitaii20ld) r 

In order to fit the mid and far-IR dust emission in the 
Grab Nebula, we constructed a grid of grain size distri- 
butions with the power-law index a ranging from 0.0-4.0, 
and Umax ranging from 0.03 to 5.0 /xm. By allowing amax 
to vary, we hope to determine to what size SN grains 
grow, and to compare the best-fit parameters to theoret- 
ical models for grain growth in Type IIP SNe. Fitting the 
minimum grain size cut-off amin allows us to estimate the 
amount of future grain processes in SNRs. In the case 
of the Grab Nebula, we fixed amin to 0.001 fim since the 
minimum grain size cut-off cannot be constrained by our 
mid-IR spectrum. In addition, since the dust in the Grab 
is radiatively heated by the PWN and the reverse shock 
has not yet reached the PWN, we do not expect that any 
grain destruction has yet occurred. Fixing amin will not 
affect our final dust mass estimate since the mass of these 
smaller grains is negli gible compared to the grains emit- 
ting i n the far-IR (see lTemim et al]|2012bl : iGomez et al.l 
[2012 ). 

2.4. Fitting of the IR Spectrum 

We fitted our entire grid of models to the observed 
dust emission from the Grab Nebula that include the av- 
erage Spitzer IRS spectrum of the dust emission scaled 
to the t otal synchrotron and line-subtracted MIPS 24 
fiui fiux (|Temim et al.ll20"l2b[ ). the synchrotron and line- 
subtracted Herschel PAGS 70, 100, and 160 fim fiux mea- 
surements and the SPIRE 250 and 350 /Ltm measurements 
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Fig. 5. — Top : Best -fit two-temperature carbon (BE) fit as used 
in lGomez et al.l 1)20121 '). The blue curve represents the warm grains 
with a temperature of 63 K, while the red curve represents the cold 
34 K grains. The black curve is the sum of the individual compo- 
nents. Middle: Our best-fit dust heating model for carbon (AC) 
grains. The rainbow colored curves represent individual spectra of 
dust grains of different sizes that are heated to a range of temper- 
atures from 35-68 K. The black curve represents the sum of the 
individual rainbow-colored spectra. Bottom: The same model as 
in the middle plot, but with the individual spectra added cumula- 
tively. 



TABLE 1 
Best-Fit Model Parameters 



Composition 


r (pc) 


a 


^max 






Silicates 


0.5 


3.5±0.1 


5.0 


f ... 


27.1 


Carbon (AC) 


0.5 


3.5±0.1 


0.6±0.2 


11.7 


Carbon (BE) 


0.5 


3.5±0.1 


> 


1.0 


30.7 



Note. — Best-fit model parameters where r is the 
distance of the dust grains from the center of the PWN, 
a is the power-law index on the grain size distribution, 
and amax is the maximum grain size cut-off. The uncer- 
tainties on the parameters for each grain composition 
were fo und from Monte Carlo simulations described in 
Section [231 

(jGomez et al.|[20T^ . where cold dust emission still con- 
tributes. The data and the best-fit model for each grain 
composition are shown in Figure [2j The best-fit distance 
r for all compositions is 0.5 pc, while the best-fit values 
for the power-law size distribution index a and the max- 
imum grain size cut-off a^ax are summarized in Table [TJ 
The dust model for the amorphous carbon (AC) grain 
composition provides the best fit to the data, since it si- 
multaneously provides a good fit to the mid-IR spectrum 
and the Herschel data points at far-IR wavelengths. The 
best-fit silicate grain composition falls somewhat short at 
far-IR wavelengths, while the carbon (BE) composition 
slightly underestimates the IRS spectrum in the mid-IR. 

2.5. Monte Carlo Simulations 

We note that the large uncertainties in the data points 
in the IRS spectrum are dominated by systematic uncer- 
tainties introduced by the subtraction of the u nderlying 
synchrotron spectrum (see lTemim et al.ll2012bD . For this 
reason, the reduced for our best fits is too small, on 
the order of ~ 0.06. Since our uncertainties are not ran- 
dom statistical uncertainties, we cannot use statistics 
in determining confidence levels for our fitted parame- 
ters. Instead, we used a Monte Carlo simulation to esti- 
mate the spread in our fitted parameters. The absolute 
values from our fits are used only to demonstrate the 
relative goodness of fit for the various models. 

In order to estimate the uncertainties in the fitted pa- 
rameters while correctly accounting for the statistical un- 
certainties in the dust spectrum introduced by the sub- 
traction of the synchrotron spectrum, we generated 5,000 
simulated spectra using the Monte Carlo method. The 
spectra were simulated by adding an average global syn- 
chrotron spectrum to the IRS data , with a synchrotro n 
power-law index of 0.42, as found bv lGomez et al.l ()2012D . 
We then produced 5,000 sample dust spectra by re- 
subtracting a global synchrotron spectrum while ran- 
domly selecting spectral indices within 3a (±0.02) of 
the best fit value of 0.42. We also randomly selected 
values for the 3.6 /xm flux from which the synchrotron 
spectrum is extrapol ated in the range of 12.65±0.63 Jy 
(|Temim et al.ll2012bl ). and we accounted from the ~5% 
uncertainty that results from normalizing the IRS dust 
spectrum to the total flux from dust in the MIPS 24 fim 
band. The Spitzer MIPS and Herschel data were simi- 
larly varied within the uncertainties listed in Table 4 of 
iGomez et al.1 ()2012D . The resulting range of simulated 
spectra are shown as the gray band in Figure |3l with the 
actual data overlaid in black. We then fltted our grid of 
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Fig. 6. — contour plots for the maximum grain size cutoff amax and the index on the power-law grain size distribution a. The contours 
indicate that an a of 3.5 and a larger grain size cutoff (> 0.1 /im) tend to produce better fits for both silicate and carbon (AC) grain 
compositions. 

dust distribution models to each of the simulated spec- 
tra to obtain a distribution in the best-fit values for the 
dust distribution power-law index a, a„iax: and the to- 
tal dust mass for each composition. The results for each 
parameter will be discussed in the following sections. 



3. RESULTS AND DISCUSSION 

3.1. Distribution of Grain Temperatures 

The main difference in using a physical dust heating 
model that allows for a distribution of grain sizes, rather 
than using two-component fits for two distinct temper- 
atures, is the resulting continuous distribution of dust 
grain temperatures that will determine the shape of the 
total dust emission spectrum. Figure 0] shows the dust 
temperature (T) as a function of grain size, assuming 
that the heating source (i.e. the Crab's PWN) is located 
at the best-fit distance of 0.5 pc. The colored curves show 
the temperature profiles for each of the three grain com- 
positions that were used in the fit. The shape of the tem- 
perature profile depends on the product L^Qabs(y,a), 
integrated of the available frequency range. An impor- 
tant thing to notice is the significant difference in tem- 
perature of the smaller grains between the carbon (AC) 
and carbon (BE) compositions. While the carbon (AC) 
grains with sizes below '-^O.Ol /xm are heated to a temper- 
ature of ^^68 K, the same-sized carbon (BE) grains are 
only heated to ~ 55 K. The large difference in dust tem- 
peratures between these two carbon compositions is not 
due to the differences in grain properties or the shape 
of the absorption coefficients Qabs, but is instead pri- 
marily due to the available wavelength coverage of Qabs- 
In Section 12. 1[ we mentioned that we only consider Li, 
up to an energy of about 0.6 keV since the energy de- 
posited in the dust above this energy becomes negligible. 
However, the absorption coefficients for the carbon (BE) 
grains are only available down to a wavelength 0.05 /im, 
which means that was only integrated up to an en- 
ergy of 0.02 keV. Accounting for energies above 0.02 keV 
is especially important for heating by PWNe, since their 
non-thermal spectra peak at these energies. The car- 
bon (AC) composition therefore provides a more phys- 
ical dust heating model for the Crab. The additional 
coverage between 0.002 and 0.05 /xm (0.02-0.6 keV) for 



the carbon (AC) grains raises the temperature for all 
grain sizes, but has a most apparent affect on the smaller 
grains (see Figure S]) . This significant difference in dust 
temperature emphasizes the importance of using grain 
compositions with adequate wavelength coverage for the 
absorption coefficient when constructing models to de- 
scribe dust emission in SNRs, especially dust surround- 
ing PWNe. 

Figure [5] shows the comparison between a two- 
temperature fit to the dust emission in the Crab Neb- 
ula using a carbon (BE) composition (top panel), as in 
iGomez et al.l (|2012[ ). and our best-fit model for the car- 
bon (AC) grains with parameters listed in Table [T] (mid- 
dle and bottom panels). The middle panel of Figure [5] 
shows spectra of grains with different sizes as individual 
rainbow colored curves that all add up to the best-fit 
black curve. The bottom panel shows the same model, 
but with every tenth cumulative sum of individual spec- 
tra shown in rainbow colors. While both models provide 
equally good fits to the IR data, the physical dust heat- 
ing model allows us to constrain physical properties of 
the dust, and significantly affects the total dust mass 
responsible for the IR emission (see Section [531) . 

3.2. Grain Size Distribution 

The dust heating model also allows us to place some 
constrains on the shape of the grain size distribution and 
the maximum size of grains that formed in the Crab Neb- 
ula. Figure |5] shows the contour plots for a, the index 
of the power-law size distribution, versus the maximum 
grain size cut-off amax- The parameter values for mod- 
els with minimum values are listed in Table [T] The 
contour plots show the relative goodness of fit for the 
range of parameter values used in our grid of dust dis- 
tribution models. For all compositions, the best-fit value 
for the power-law index is a=3.5. The fits to the 5,000 
Monte Carlo simulations show that the power-law index 
is very well constrained with a standard deviation of only 
0.1. This is t he same size distrib ution that is found for 
the ISM (e.g. IMathis et al.|[T977l) . Physical modeling of 
dust emission for a larger sample of SNRs is needed to 
determine if similar size distributions are found in other 
core-colla pse SNe. 

lu iTemim et al.l (|2012b[ ). we found that the dust emis- 
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TABLE 2 

Total Dust Mass (M©) 



Composition 



Observations Nucleosynthesis Models 

This Work Gomez et al. (2012) WW95 N07 WH07 



Silicates 0.128 ±0.004 

Carbon (AC) 0.027 ± 0.003 
Carbon (BE) 0.050 ± 0.002 



24+°-32 
0.11 ±0.01 



0.08 
0.05 
0.05 



0.32 
0.10 
0.10 



0.12 
0.11 
0.11 



Note. — . The uncertainties on the dust mass from this work were found 
from Monte Carlo simulations described in Section 12.51 The dust mass listed 
for the nucleosynthesis models is the maximum allowed dust mass assuming 
a 100% co ndensation efficiency, and the yields for a 11 Mq p rogenitor from 
IWooslev"^ Weaven U995|') (WW95). and a 13 Mq progenitor from lNomoto et aP 
II2006I) fN07) and lWooslev fc Heeej ||2007| ) (WH07). 



sion in the Crab observed by Spitzer implies a small grain 
size of < 0.05 /xm. However, since we only had cov- 
erage up to 70 fim, we were sampling only the smaller 
grains that were heated to temperatures of around 50 K 
(see Figure HJ . The addition of the Herschel long wave- 
length data has now allowed us to place constraints on 
the maximum size of dust grains in the Crab's filaments. 
Our best-fit model appears to favor a fairly large max- 
imum grain size cut-off (see Figure [B]) ■ The value for 
o-max for the best-fit model is 0.6 fim or larger, depend- 
ing on the grain composition (Table [Ij . This is actu- 
ally more consi stent with models of dust formation i n 
Type IIP SNe feozasa et al.|[2009t iNoiawa et al.|[20Tol) . 
These models show that the dust mass of grains formed 
in Type IIP SNe with massive H-envelopes is dominated 
by grains larger than 0.03 /im, while the mass of dust 
for Type lib SNe is dominated by very small grain < 
0.006 fim. Probing the dust grain size in SNRs can 
therefore be an extremely useful tool for investigating 
the properties of the SN progenitor. The Crab Nebula is 
thought to be a result of a Type IIP explosion, primar- 
ily based on the low expansion velocities of the ejecta 
filaments and the low progenitor mass (e.g. IChevalieil 
^005: M acAlpine fc Sattc rficld 2008). The large dust 
grain radii inferred from our models are consistent with a 
Type IIP origin. We note that the larger grain size also 
implies that a large fraction of the dust m ass will sur- 
vive the even t ual reverse shock in teraction ("Dwe kl 120051: 
Kozasa et al.l 120091 : IBianchi et al., i2009> .Nozawa et al.l 
20101) . 



3.3. Revised Dust Mass for the Crab Nebula 

The total dust mass for the Crab Nebula that re- 
sults from our heating model is 0.128 ± 0.004 Mq and 
0.027 ± 0.003 M0 for silicate and carbon (AC) grains, re- 
spectively (Tabled]). Figure [7] shows the histograms of 
the total dust mass for the best-fit model for each of the 
5,000 simulated spectra. A Gaussian distribution was 
fitted to the histogram and its mean and a taken as the 
estimated mass and its uncertainty. The plots shows that 
the total mass is very well constrained for both silicate 
and carbon grain compositions. We note that the total 
mass for each composition is significant lower than the 
masses estimated from the two-component temperature 
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+0.32 
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Mq for silicates and 0.11 ± 0.01 Mq for 



o 



carbon (|Gomez et al.ll2012[ ) (see Tabled. The silicate 
and carbon (BE) compositions that we used in ou r fits 
are the same as those used in (I Gomez eFaLllMl . and 
we confirmed this by reproducing their total dust masses 



1 500 \ 



o 1000 



Silicates ■ 



500 




.12 0.15 0.-4 

Total Dust Mass (Mg) 



0.15 



?000 [ 



13 1 500 



1000 - 



Carbon (ac) . 




500 



0.015 0.020 0.025 0.050 0.055 

Total Dust Mass (Mg) 



0.040 



Fig. 7. — Distribution of total dust masses from the fits to 5000 
Monte Carlo simulations of the IR data using silicate (top) and 
carbon (AC) grains (bottom). The best-fit Gaussian is overplotted 
on the histograms as a solid black line. We use the mean and 
standard diviation of the Gussian for the best-fit total mass and 
uncertainty, 0.128 ± 0.004 Mq for silicates and 0.027 ± 0.003 Mq 
for carbon grains. 

with a two-temperature dust fit. We find that if we use 
these same absorption coefficients in our dust heating 
model, the total mass of dust is reduced by a factor of 
r^2. The total mass for carbon grains is lowered even 
further by using more realistic Qabs values with cover- 
age at lower wavelengths. As discussed in Section 13.11 
using the carbon (AC) grains raises the temperature of 
the smaller grains, which will in turn contribute more 
flux to the spectrum and lower the required number of 
large cooler grains that dominate the mass. The signifi- 
cantly different mass estimates are therefore the results 
of a combination of a physical grain size and temperature 
distribution, and a selection of well measured absorption 
coefficients that will account for the absorption of higher 
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energy photons. 

The fact that carbon (AC) grains provide a signifi- 
cantly better fit to the dust emission in the Crab Nebula 
suggests that the dust mass may either be dominated 
by carbon grains, or possibly a mixture of carbon and 
silicate grains. The total mass estimated from silicates 
should therefore be treated as an upper limit, and the 
lower carbon mass of 0.027 ± 0.003 as a more realistic 
estimate for the total dust mass in the Crab Nebula. This 
lower dust mass is further supported by nuclcosynthetic 
yields from core- collapse SNe for progenitors in the 11-13 
r ange (iWoosley fc Weaveri 119951: iWooslev fc Hegerl 
I2007t iNomoto et al.ll2006D . We assumed a 100% grain 
condensation efficiency and calculated the upper limit 
on the carbon and silicate grain mass from the total 
yield of carbon, oxygen, magnesium, silicon, and iron. 
To obtain the upper limit on the silicate grain mass, we 
summed the maximum possible masses for Si02, MgO, 
and Fe304 grains. The results ar e listed in Tabled and 
they show that the iGomez et al.l (|2012f ) mass estimates 
are clearly too large, especially considering that the pro- 
genitor of t he Crab Nebula is e stimated to be around 
- 9.5 Mq (IMacAlpine fc Satter ficld 200^. The total 
mass of silicate dust derived from our models is also too 
large compared to nucleosynthesis yields, suggesting that 
the dust in the Crab Nebula is mostly in the form of car- 
bon grains. This work implies that a physical dust model 
that allows for continues size and temperature distribu- 
tions will revise the total dust mass estimates for many 
SNRs, especially the recent estimates made from fitting 
simplified two-component models to account for the far- 
IR emission, including the contr oversially high mass o f 
0.4-0. 7Mq found for SN 1987 A (IMatsuura et al.l[2QTTh . 

4. CONCLUSIONS 

We modeled the mid and far-IR dust emission from the 
Crab Nebula in order to constrain the dust grain prop- 
erties and total dust mass. We use a physical dust heat- 
ing model in which a continuous power-law size distri- 
bution of dust grains is radiatively heated by the PWN. 



The best-fit models favors a grain size distribution with 
a power-law index of 3.5, and a relatively large maxi- 
mum grain size cut-off of > 0.6 //m, consistent with the- 
oretical predictions of dust formation in Type IIP SNe 
(jKozasa et al.l [20091: iNoiawa et al.llMol) . The IR emis- 
sion is best described by an amorphous carbon grain 
composition with a total mass of 0.027 ± 0.003 Mq, sig- 
nificantly lower than previous estimates from more sim- 
plified two-temperature models. The difference in the 
dust mass is due to a combination of fitting a continu- 
ous grain size and temperature distribution, and apply- 
ing grain absorption coefficients with a more complete 
wavelength coverage that allows us to account for heat- 
ing from higher energy photons from the PWN. 

In recent years, most large dust mass estimates for 
SNRs assumed the simplified one or two-temperature 
models, including the mass estimates from Herschel 
data for Cas A, SN 1 987 A , and t he Crab Neb- 
ula fe.g. iBarlow et al.l [2OIOI: Sibthorpe et al.l I20T0I: 
IMatsuura et al.ll2011l:lGomez et al.i,2012.) . Our modeling 
results imply that a more realistic dust emission model 
will revise, and may significantly reduce the mass esti- 
mates for many SNRs. The amount of SN-condensed 
dust in the Crab Nebula is significantly lower than that 
required to produce the massive amount of dust observed 
in high-reds hift galaxies even if the dust is not destroyed 
in the ISM ()Dwek et al.ll2007[ ). 

This will have important implications on the ques- 
tions of dust production in SNe and its ultimate survival. 
While the model presented in this paper assumes radia- 
tively heated dust, we are also combining varying grain 
size distributions with a dust heating model for collision- 
ally he ated grains (see lArendt et al.l 120101 : iTemim et al.l 
I2012al) . We plan to apply our models to revise dust mass 
estimates for a large number of PWNe and SNRs, and 
use them to constrain dust grain properties that may 
shed light on the nature of the SN progenitors. 

We thank Rick Arendt and George Sonneborn for use- 
ful discussions and comments on the paper. 
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